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ABSTRACT 

Several glitches have been observed in young, isolated radio pulsars, while a clear detection in accretion-powered X-ray pulsars is 
still lacking. We use the Pizzochero snowplow model for pulsar glitches as well as starquake models to determine for the first time 
the expected properties of glitches in accreting pulsars and their observability. Since some accreting pulsars show accretion-induced 
long-term spin-up, we also investigate the possibility that anti-glitches occur in these stars. We find that glitches caused by quakes in 
a slow accreting neutron star are very rare and their detection extremely unlikely. On the contrary, glitches and anti-glitches caused 
by a transfer of angular momentum between the superfluid neutron vortices and the non-superfluid component may take place in 
accreting pulsars more often. We calculate the maximum jump in angular velocity of an anti-glitch and we find that it is expected to 
be MV g , » 1(T 5 - 1(T 4 rad s _1 . We also note that since accreting pulsars usually have rotational angular velocities lower than those of 
isolated glitching pulsars, both glitches and anti-glitches are expected to have long rise and recovery timescales compared to isolated 
glitching pulsars, with glitches and anti-glitches appearing as a simple step in angular velocity. Among accreting pulsars, we find that 
GX 1+4 is the best candidate for the detection of glitches with currently operating X-ray instruments and future missions such as the 
proposed Large Observatory for X-ray Timing (LOFT). 

Key words, dense matter - pulsars: general - stars: neutron - X-ray: binaries - pulsars: individual: GX 1+4 


1. Introduction 


Isolated pulsars are considered remarkably stable rotators show¬ 
ing long-term spin-down caused by the emission of electromag¬ 
netic dipole radiation, relativistic outflows, and possibly gravita¬ 
tional waves (Aasi et al. 2015). However, they have sometimes 
been observed to glitch. Glitches are sudden increases in the 
spin rate followed by a relaxation toward a steadily, long-term, 
decreasing spin rate. They have been observed in over 100 iso¬ 
lated radio pulsars and magnetars (see, e.g., Espinoza et al. 2011; 
Yu et al. 2013; Dib et al. 2008), with jumps in angular velocity 
up to AQ g i « 10~ 4 rad s -1 . 


There are two main models to explain glitches: the star- 
quake glitch model (e.g., Baym & Pines 1971) and the super¬ 
fluid vortex unpinning model (Anderson & Itoh 1975; see also 
Haskell & Melatos 2015 for a recent review on models of pulsar 
glitches). In the starquake glitch model of Baym & Pines (1971), 
a rapidly rotating pulsar has an oblate crust that deforms toward 
an almost spherical shape as the pulsar slows down. This leads 
to a sudden crack in the crust and a decrease in the moment of 
inertia resulting in a sudden increase in the angular velocity of 
the star. In the model proposed by Anderson & Itoh (1975) the 
glitch is caused by a sudden unpinning of neutron superfluid vor¬ 
tices from the lattice nuclei in different regions of the interior of 
a neutron star (NS). The angular momentum stored in the super¬ 
fluid vortices is then transferred to the non-superfluid component 


of the star, leading to the observed fast spin-up of the NS rotation 
angular velocity. 

No clear observational evidence that accreting NSs experi¬ 
ence glitches has been presented so far. Nonetheless, they have 
been proposed as a possible explanation of the peculiar timing 
behavior of some accreting pulsars. Galloway et al. (2004) ob¬ 
served a sudden spin-up (with a fractional change in frequency 
of ~ 3.7 x 10 5 in < 10 hr) in the accreting 18.7 s pulsar 
KS 1947+300 with RXTE, that was interpreted as a glitch event. 
Contrary to what happens in the classical glitch scenario, the 
sudden spin-up was observed when the NS was already spin¬ 
ning up because of the transfer of angular momentum from the 
accretion disk. Galloway et al. (2004) did not rule out the pos¬ 
sibility that the sudden spin-up was caused by a short episode 
of enhanced accretion. Klochkov et al. (2009) studied the pulse 
period variations of Her X-1 using SwiftfB AT data. They found 
two spin-down episodes followed by spin-up intervals where the 
spin period increased exponentially, a phenomenon typical of the 
post-glitch recovery stage, but with an opposite sign. 

Observations of glitching pulsars showed that the glitch 
frequency increases with the long-term spin-down rate 
(Espinoza et al. 2011). Since the spin-down caused by elec¬ 
tromagnetic braking decreases as pulsars age, the detection of 
glitches in old pulsars (> 10 4 yr, Espinoza et al. 2011), such as 
accreting pulsars in X-ray binaries, is expected to drop. In this 
paper we show that glitch frequency in accreting pulsars in X- 
ray binaries can actually be higher than expected for old pulsars. 
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Table 1 . Rotational parameters and luminosities of four accreting pulsars showing among the largest and stable spin-down rates (f2„) in XRBs. £1 
is the rotational angular velocity, r„ is the torque reversal timescale, and L x is the X-ray luminosity. 


Name 

G 

rad s~ 

1 

Goo 

rad s~ 1 2 


T OO 

d 


L x 

erg s _1 


GX 1+4 

~ 0.042 

(1) 

-3 x 10- 11 

(1) 

> 9900 

(1) 

IO 35 - 10 36 

(1) 

OAO 1657-415 

* 0.17 

(2) 

-1.3 x 10’ 11 

(2) 

* 200 

(3) 

10 36 - 10 37 

(2) 

4U 1626-67 

* 0.8 

(4) 

-4.4 x 10^ 12 

(4) 

« 6600 

(5) 

> io 36 

(4) 

4U 1907+09 

« 0.014 

(6) 

-2.3 x 10~ 13 

(6) 

> 5500 

(6) 

10 36 

(7) 


References: 

(1) Gonzalez-Galan et al. (2012); (2) Jenke et al. (2012); (3) Barnstedt et al. (2008); (4) Chakrabarty et al. (1997); (5) Beri et al. 
(2014); (6) Inam et al. (2009); (7) in’t Zand et al. (1997). 


The aim of this paper is to assess the conditions under which 
glitches are more likely to occur in an accreting pulsar and in¬ 
vestigate, for the first time, the expected properties of glitches in 
these objects. 

In Sect. 2 we calculate the expected jump in angular velocity 
and the interglitch interval in the starquake and superfluid vortex 
scenarios, assuming that glitches occur during a long-term spin- 
down. Then, we investigate the properties of anti-glitches (sud¬ 
den spin-down) during long-term spin-up of the NS. We mod¬ 
ified the model of Pizzochero (2011) to calculate the jump in 
angular velocity in this scenario. In Sect. 3 we discuss the ob¬ 
servability of glitches and anti-glitches in accreting X-ray pul¬ 
sars. 

2. Properties of glitches in accreting pulsars 

Pulsars in X-ray binary systems (XRBs) can experience spin- 
up and spin-down caused by the interaction between the ac¬ 
cretion flow and the magnetosphere of the NS (Pringle & Rees 
1972; Rappaport & Joss 1977). The interaction may occur di¬ 
rectly from the wind of the donor star or through an accretion 
disk (e.g., Bildsten et al. 1997 and references therein). In some 
XRBs (see examples in Table 1) the accretion torque is responsi¬ 
ble for variations of the angular velocity of the NS of magnitude 
comparable to that caused by the electromagnetic braking torque 
in young, glitching pulsars (10~ 15 - IO -10 rad s -2 , Lyne et al. 
2000 and, for comparison. Table 1). Since the rate of glitches in¬ 
creases with the spin-down rate (see following sections), glitches 
in XRBs should be quite frequent. In the next sections we derive 
their main properties in different scenarios. 

2.1. Starquake glitch scenario 

The Baym-Pines model allows the time Af q to the next quake to 
be estimated from the magnitude of the preceeding one through 
the equation 

2A 2 

^ = r ^ |Af| * (1) 

where T = -G/EG is the characteristic timescale at which the 
pulsar slows down because of the loss of rotational energy, G 
is the rotational angular velocity of the NS, is the long-term 
spin-down rate, | Ae| is the reduction in oblateness e caused by the 
quake in the previous glitch, A and B are coefficients describing 
respectively the gravitational and elastic energy stored in the NS 
as a result of its rotation, and Iq is the moment of inertia. The 


values of A, B, and Io depend on the equation of state (EoS) 
of the neutron star core and on the crust model adopted. Here 
we adopt the parameters obtained by Pandharipande et al. (1976) 
and Zdunik et al. (2008), which updated the Baym-Pines model 
using different EoSs and crust models. 

We assume different values of Ae ranging from 10 10 to 10 6 
(which are the typical values of glitches observed in isolated pul¬ 
sars). We use G and fl„ of the four XRBs in Table 1 . We find that 
for any value of G, EG, A, B, Iq, and Ae, A t q is always greater 
than 10 s yr. The value of Af q can be <sc 10 s yr for Ae <s 10 l0 . 
However, starquakes with such small Ae would be difficult to 
detect. Therefore, the detection of a glitch produced by a star¬ 
quake in an accreting pulsar is extremely unlikely. Moreover, 
Zdunik et al. (2008) showed that an accreted crust is softer to 
an elastic deformation than the crust of an isolated pulsar 1 . As a 
consequence, Af q would be longer and glitches caused by quakes 
rarer. 

2.2. Superfluid vortex scenario 

Our calculation of the interval between glitches Af g i and of the 
jump in angular velocity AG g i of accreting pulsars in the su¬ 
perfluid vortex scenario is based on the “snowplow” model of 
Pizzochero (2011). In this model, the matter of the NS is di¬ 
vided into two components: the neutron superfluid and the nor¬ 
mal component (charged components that corotate with the pul¬ 
sar magnetic field). In the core (r < R c « 9.3 km) and inner crust 
( R c < r < R — 10 km) the rotating superfluid is organized as an 
array of vortices parallel to the spin axis of the NS 2 . The vortices 
are pinned to the crustal lattice of ions, but are not coupled with 
the normal component of the star. Therefore, although the crust 
spins-down because of the electromagnetic braking, the super- 
fluid conserves its angular momentum 3 * * . As the NS spins down, 
a rotational lag builds up between the superfluid vortices and the 
normal component. When it reaches a critical value, a hydrody- 
namical force acting on vortices, the Magnus force, unpins and 
moves them out. The Pizzochero (2011) model introduces for the 

1 Zdunik et al. (2008) pointed out that the calculations for the elastic 
properties of accreted crusts are subject to high uncertainties, and there¬ 
fore they must be taken with caution. 

2 The value R c = 9.3 km is obtained assuming a n = 1 polytropic 
density profile to describe the core and inner crust, with R = 10 km and 
M = 1.4 M 0 for the radius of the inner crust and mass of the neutron 
star (see Pizzochero 2011). 

3 The angular velocity of the superfluid component G s is proportional 

to the number of vortices. As long as the vortices remain pinned to the 

lattice of ions, their number is conserved, hence fl s also does not vary. 
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first time a density profile of the pinning force with maximum 
value f m = 1.1 x 10 15 dyn cm -1 at densities p m ~ 0.2po, where 
po = 2.8 x 10 14 g cm ~ 3 is the nuclear saturation density. Because 
of the shape of the density profile of the pinning force, the criti¬ 
cal rotational lag for depinning is maximum at x m ~ 9.7 km (in 
cylindrical coordinates, with the z-axis parallel to the spin axis). 
Therefore, vortices from x < x m accumulates in a vortex sheet at 
x m . When the rotational lag at x m reaches the critical value, the 
vortex sheet suddenly moves out and exchanges the stored angu¬ 
lar momentum with the normal component, causing a glitch. 

With this model it is possible to determine the fraction of vor¬ 
tices pinned to the lattice nuclei in different regions of the inte¬ 
rior of a neutron star and the timescale of the pinning-unpinning 
process. The model can be used to predict the interval between 
glitches Afgi, the jump in angular velocity AfL] and the jump in 
angular acceleration AO g i during a glitch. 

We leave the observables il and Ooo as well as the two free 
parameters of the model f m and Y s \ free to vary. The parameter 
describes the fraction of vorticity coupled to the normal crust 
on timescales of the glitch rise time (< 40 s; Dodson et al. 2002). 

According to Pizzochero (2011), the interval between 
glitches can be expressed as 


Afgi — 28.8 


R l 


/l5 


M\a |Q_ n 


yr, 


( 2 ) 


where R 6 = R NS /10 6 * * cm, Mia = M NS /1.4M 0 , f] 5 = 
/m/10 15 dyn cm -1 , and 0 _n = Goo/ICT 11 rad s -2 . 

When Q m increases, At .a decreases resulting in a larger prob¬ 
ability that a glitch will be observed. Indeed, the greater the spin- 
down rate of the crust, the more quickly the critical rotational 
lag for depinning will be reached. In addition, Eq. (2) shows that 
when f m increases. At a also increases owing to the larger un¬ 
pinning critical lag required by the Magnus forces acting on the 
neutron vortices to overcome the pinning force. 

The jump in angular velocity is given by (Pizzochero 201 1), 


AQgi =* 1.25 x 10“ 


2o.95^6/l5 


M\a[\ - (2o.95(l _ I 0 . 05 )] 


■ rad s 


(3) 


where 60.95 = 6/0.95 (Zuo et al. 2004), Q is the superfluid frac¬ 
tion of the star, and T 0.05 = Tgi/0.05 (Pizzochero 201 1). 

When / m increases, AQ gl increases because of the larger 
number of neutron vortices unpinned simultaneously (see 
Pizzochero 2011). Equation (3) also shows that when Y„\ in¬ 
creases, Afi g i decreases because of the reduced number of vor¬ 
tices storing angular momentum coupled to the normal compo¬ 
nent of the neutron star. 

Since accreting pulsars usually have angular velocities 
smaller than those of isolated pulsars, we expect the glitch rise 
time and the recovery time after glitch 4 * to be longer in accret¬ 
ing pulsars. Indeed, these timescales depend on the coupling 
timescale between the superfluid and the normal component, 
which is inversely proportional to the angular velocity of the NS: 


T OC 


1 

Q 


(4) 


(see, e.g., Haskell & Antonopoulou 2014 and references 
therein). Therefore, since the glitch rise time of Vela pulsar 

4 The post-glitch recovery stage, which follows the sudden spin-up, 

is characterized by a quasi-exponential relaxation of the angular ve¬ 

locity time derivative of the star toward the spin-down rate which the 

pulsar would have had in the absence of the glitch (Baym et al. 1969; 

Alpar et al. 1984; Jones 1990). 



Fig. 1. Schematic representation of the cylindrical geometry assumed 
in Sect. 2.3. The z-axis is parallel to the spin-axis of the NS. The dotted 
line shows a vortex line with half length sjR 2 — x 2 . 


(£2vela = 70 rad s -1 ) is rveta ~ 1 - 10 s (Haskell et al. 2012), 
in a slow rotator such as GX 1+4 it would be Text +4 ~ 10 3 - 
10 4 s. Similarly, the recovery timescales of slow rotators become 
longer as well. The glitch is expected to appear as a simple step 
in angular velocity. 

Glitches with relaxation timescales longer than those of Vela 
pulsar have been observed in some slower anomalous X-ray pul¬ 
sars (Haskell & Antonopoulou 2014; Dib et al. 2008). The glitch 
rise time has not yet been measured. An upper-limit of ~ 40 s has 
been placed by Dodson et al. (2002) for Vela pulsar. 

2.3. Anti-glitch scenario 

Some XRBs show long-term spin-ups with timescales on the or¬ 
der of ~ 10 years (e.g., 12yrof4U 1626-67; Chakrabarty et al. 
1997 and referencese therein), and rates up to Cl « 10 10 rad s 2 
(e.g., KS 1947+300; Galloway et al. 2004). In principle they are 
good candidates for experiencing the anti-glitch scenario pro¬ 
posed by Pines et al. (1980), in which a pulsar can show a sud¬ 
den spin-down caused by a mechanism of angular momentum 
transfer similar to that of glitches. 

Here we adapt the snowplow model of Pizzochero (2011) to 
the scenario proposed by Pines et al. (1980) to calculate the ex¬ 
pected jump in angular velocity caused by an anti-glitch AQ a „|. 
Hereafter we will use cylindrical coordinates as shown in Fig. 1. 

In the glitch scenario, when the pulsar slows down, vortices 
from the core region are expelled outward, allowing the angular 
momentum of the core to rebalance easily, while the vortices in 
excess are accumulated at the boundary x c and then expelled out¬ 
ward as vortex sheet. In the anti-glitch scenario, on the contrary, 
the crust accelerates. Before the anti-glitch, new vortices are cre¬ 
ated at R , near the equator, as the star accelerates. They are ac- 
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cumulated in a vortex sheet by the external depinning front 5 that 
moves inward from R to x m . In the meantime, as the internal 
depinning front moves across the region x < x m , the vorticity 
moves from the inner crust to the core and rebalances the angu¬ 
lar momentum in the core region. Since the number of vortices 
within x m does not change, their accumulation in the core region 
must correspond to a depletion in the regions of the inner crust 
with lower pinning potential, which is around x c , the boundary 
between the two regions (see Pizzochero 2011 and figures 1 and 
3 therein). We assume that when the vortices accumulated at x m 
by the external depinning front are released simultaneously, they 
fill the depleted region around x c . Hence, the transfer of angular 
momentum will take place in the region x c < x < x m . We point 
out that the minimum (x m ; n ) and maximum (x max ) radii of the 
depletion region are unknown. A precise estimate of these val¬ 
ues requires a detailed study of vortex dynamics. Nonetheless, 
it is still possible to roughly estimate the order of magnitude of 
the jump in angular velocity caused by an anti-glitch relying on a 
few reasonable assumptions on x^ and x max as described below. 

According to Seveso et al. (2012), the angular velocity of the 
superfluid component of the star £2 s (x) is given by 


O s (x) 


k N(x) 
2n x 2 


(5) 


where N(x) is the number of vortices in a cylindrical region of 
radius x, k — nh/m^, and ;«n is the mass of a neutron. Using Eq. 
(5), the number of vortices removed from the region x m < x < R 
and accumulated at x m before the anti-glitch is 


A v = — (R 2 - x^)AQ cr (x m ), 


(6) 


where AQ cr (x m ) is the critical rotational lag for depinning at x m . 
Since we do not know the exact size of the re-coupling region, 
we calculate the angular momentum released during the anti¬ 
glitch AL a - g i based on two different assumptions. In the first 
case we assume that N v vortices will be uniformly distributed in 
x min < x < x m by integrating N(x) = A v (x 2 - x 2 jn )/(x^ - x 2 J, 
where x m ; n corresponds to the minimum radius of the depletion 
zone around the boundary x c . From the definition dL = O s t// s , 
AL a _ g i can be expressed as 


r Xm x 2 — x 2 ■ f"JR 2 --c 2 

AL a _ gl = 2 kN v - - ^xdx p s (r)dz , (7) 

Jx mm *m - Jo 


5 At any time t before the glitch, the external and internal depinning 
fronts represent the radial distance x ext (r) and Xi„ t (f) (with x ext (0 > x m 
for the external depinning front and x; nt (t) < x m for the internal depin¬ 
ning front) where the rotational lag AQ = O s (x) - Q between the local 
superfluid angular velocity fl s (x) and that of the normal component Q 
is equal to the critical lag for depinning Afl cr . As the star spins up, a 
rotational lag builds up between the superfluid vortices and the normal 
component. The critical lag for depinning is highest at x m . Therefore, 
the external depinning front moves inward and the internal depinning 
front outward, both to x m . Since the NS spins up, the relative velocity 
vl - v s of the vortex lines pinned to the lattice (vl) with respect to the 
superfluid (v s ) increases. The Magnus force acting on the vortices is 
given by 

F m = kpA x (v L - v s ), 

where e z is the unit vector pointing along the rotation axis. Because of 
the direction of F m , the vortices unpinned by the external and internal 
depinning fronts move inward (while in the glitch scenario, they move 
outward). 


where r = Vx 2 + z 2 and 


Ps(r) = Q 


nM sin(nr/R ) 
W n r/R 


is the superfluid density profile (Pizzochero 2011). In the second 
case, we assume that N v vortices will be uniformly distributed in 
x c < x < x max (x max < x m ), where the pinning potential is lower. 
In this case, AL a _ g j is 


r dSr-X 1 

N(x)xdx J p s (r)dz , 


AL a _ g i - 
where 

N(x) =N v x 


(x 2 - xj?)/(x; 

1 


max 


X 2 ) 


if x c A x A x max 
if x max < x < x m 


( 8 ) 

(9) 


The jump in angular velocity of the normal component of the NS 
caused by the anti-glitch is thus given by 


AQ a g i — 


AL a _ g i 

/,ot[l - 2d - Egl)] ’ 


( 10 ) 


where 7 tot is the total momentum of inertia of the NS. In Eqs. (7) 
and (8) we integrate on the cylindrical regions x m j n < x < x m 
and x c < x < x max , where x mm and x max are unknown. Un¬ 
der the reasonable assumption that x m j n and x max are close to 
x c , where the pinning potential of the inner crust is lower, we 
find that |AQ a _ g i| is approximately 10~ 5 - 10 4 rad s -1 assuming 
f m = 1.1 x 10 15 dyncm 1 and T g i = 0.05. 

We note that Eq. (2) can also be used to describe the interval 
between anti-glitches. Indeed, Af g i is the time required to reach 
the critical rotational lag for depinning, no matter whether the 
rotational lag between the superfluid vortices and the normal 
component has been built through a long-term spin-down or a 
long-term spin-up. Hereafter we will use Af g i both for glitches 
and for anti-glitches. 


3. Observability 

The results of the previous section indicate that, in the case of 
accreting pulsars, glitches (anti-glitches) should have maximum 
jumps in angular velocity Afl g i » 1.3 X 1CT 4 rad s -1 (AQ a _ g i w 
1(U 5 - 1(T 4 rad s _1 ), rise time Tri se ~ (10 2 - 10 3 )fi _1 s, and 
long recovery timescales making them appear as simple steps in 
angular velocity leaving O nearly unchanged. 

In accreting pulsars, accretion torque variations (including 
sign reversals) result in changes in the NS angular velocity. Sev¬ 
eral transitions between spin-up and spin-down have been ob¬ 
served (Nelson et al. 1997) and have been explained so far with 
sharp changes in the coupling between the matter flowing to¬ 
ward the NS and the magnetosphere. Several models have been 
proposed both for disk- and wind-fed pulsars (e.g., Zhang & Li 
2010 and references therein; Wang 1981). However, a widely 
accepted model does not yet exist. The timescale too between 
two torque reversals varies from source to source. Some pul¬ 
sars show numerous torque reversals with moderate spin pe¬ 
riod derivative (e.g.. Vela X-l, Her X-l; Bildsten et al. 1997), 
others show quasi-periodic torque reversals superimposed on a 
long-term spin-up trend (e.g., OAO 1657-415, Barnstedt et al. 
2008; EXO 2030+375, Fermi/ GBM 6 ; SAX J2103.5+4545, 

6 The Fermi Team provides the pulse frequency history of the 
accreting pulsars monitored with the Gamma-ray Burst Monitor 
(GBM) on board the Fermi satellite since its launch (June 11, 
2008). Preliminary results are available for download from the website 
http://gammaray.nsstc.nasa.gov/gbm/science/pulsars.html 
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Camera et al. 2014), while others show steady spin-up or 
spin-down sporadically interrupted by a torque reversal (e.g., 
GX 1+4, Gonzalez-Galan et al. 2012; 4U 1626-67, Beri et al. 
2014; 4U 1907+09, Inam et al. 2009). 

It is interesting to compare, for each pulsar, too with Af g i, the 
time needed to reach the critical rotational lag for depinning. The 
accreting pulsars with sufficiently long time intervals of spin- 
down or spin-up (i.e., with too > Af g i) are those more likely to 
experience glitches or anti-glitches. In Fig. 2 we show t m and 
Afgi for a sample of accreting pulsars. We calculate Af g i from 
Geo using Eq. (2) and obtain too and Qco from previous works 
or from the Fermi /GBM monitoring archive. Black symbols re¬ 
fer to the values observed during a long-term spin-down, while 
those obtained during a long-term spin-up are plotted in blue (for 
some pulsars, e.g., 4U 1626-67, both spin-up and spin-down 
episodes have been observed, hence both cases are indicated). 
For those pulsars that exhibit quasi-periodic torque reversals su¬ 
perimposed on a long-term spin-up, namely OAO 1657-415, 
EXO 2030+375, and SAX J2103.5+4545, we use too and Af g i 
related to the long-term spin-up. For persistent pulsars that show 
a “random walk” of the spin period with several torque rever¬ 
sals and moderate |Qco| (VelaX-1, Cen X-3, Her X-l), we give 
the lower-limits on Af„|. We show the lower-limits on t a , for pul¬ 
sars with spin-up or spin-down that are still ongoing (GX 1+4, 
OAO 1657-415, SAX J2103.5+4545, spin-up of 4U 1626-67, 
and EXO 2030+375). GX 1+4 has too > Af g i, making it the 
best candidate for observing glitches. Other good candidates are 
4U 1626-67 and OAO 1657-415, with too ~ Ar g ]. We point out 
that, in principle, all the accreting pulsars can experience glitches 
or anti-glitches, although with a lower probability. Indeed, they 
have long spin period evolutions (10 4 - 10 6 yr) during which 
they might have built up a rotational lag for depinning close to 
the critical value. 

While glitches are traditionally observed with radio obser¬ 
vations, the numerous examples of glitches seen in magnetars 
clearly show that they can also be detected by X-ray observa¬ 
tories (e.g., Kaspi et al. 2000). The size of the jumps in angular 
velocity observed in the magnetars is the same or smaller than 
that expected in accreting pulsars (Afi ~ 10 7 - 10 4 rad s _1 , 
§a§maz Mug & Gogiig 2013 and references therein). 

A suitable spacing of the observations, allowing phase- 
connected timing ephemeris to be derived, is required in order to 
detect glitches and to distinguish them from other timing irregu¬ 
larities induced by variations in the accretion torque. The obser¬ 
vation of correlated changes in the source flux should help to rec¬ 
ognize accretion-induced torque variations. Existing instruments 
with a wide field of view, such as the Fermi/ GBM, Swift/ BAT, 
and INTEGRAL/IBIS are already providing repeated observa¬ 
tions of the brightest X-ray pulsars that could be used to search 
for glitches, but the observing strategy of these satellites has not 
been optimized in this respect. In the future, a relevant contribu¬ 
tion in this field could be given by the LOFT mission, proposed 
in the context of the European Space Agency science program 
(Feroci et al. 2014). The Large Area Detector (LAD) instrument 
on LOFT will be able to detect the pulsations of bright sources 
like 4U 1626-67 or GX 1+4 in observations of only a few kilo- 
seconds. A program of suitably spaced monitoring observations 
will be easily implemented, also exploiting the frequent cover¬ 
age of the Galactic plane with the LOFT Wide Field X-ray mon¬ 
itor. 
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Fig. 2. Afgi - too diagram of 13 accreting pulsars. Black refers 
to spin-down, blue to spin-up. Pulsars with horizontal arrows ex¬ 
perience a “random walk” of the spin period, and so the value 
Afgi reported in the diagram is a lower-limit. Pulsars with verti¬ 
cal arrows have ongoing spin-up or spin-down. We obtained f«, and 
AfgilDoo) from previous works (OAO 1657-415: Jenke et al. 2012; 
4U 1626-67: Chakrabarty et al. 1997; Beri et al. 2014; 4U 1907+09: 
Inam et al. 2009; SAX J2103.5+4545: Camera et al. 2014), or from 
the Fermi /GBM monitoring archive (GX 1+4, EXO 2030+375, 
2S 1417-624, RX J0520.5-6932, GRO J1744-28, KS J1947+300, 
Cen X-3, Vela X-l, Her X-1). 


4. Conclusions 

We use starquake and superfluid vortex models to outline for 
the first time the expected observational properties of glitches in 
accreting pulsars. 

We find that glitches caused by quakes of the crust in an ac¬ 
creting slow pulsar are very rare and their detection unlikely. On 
the contrary, glitches caused by the transfer of angular momen¬ 
tum between the superfluid and the normal components may take 
place more often. In the superfluid vortex scenario they can have 
maximum jump in angular velocity AQ g i ~ 1.3 x 10 4 rad s 1 . 
Some accreting pulsars show long-term spin-up periods during 
which they may undergo anti-glitches. Therefore, we modify the 
snowplow model of Pizzochero (2011) on the basis of the anti¬ 
glitch scenario proposed by Pines et al. (1980) to determine the 
expected properties of anti-glitches in these objects. We find that 
anti-glitches in accreting pulsars have maximumjumps in angu¬ 
lar velocity of about the same order of magnitude or ten times 
smaller than AO g i. Both glitches and anti-glitches are expected 
to have long rise times (rri se ~ 1 () 2 - 10 ’fl 1 s) and long recovery 
timescales, with the glitch appearing as a simple step in angular 
velocity. 

By comparing the timescales between two consecutive 
torque reversals too and the interglitch timescale Af a i of a sam¬ 
ple of accreting pulsars, we find that GX 1+4 has too > Af g i (see 
Fig. 2). It is therefore the best candidate for observing glitches. 
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Other good candidates with too ~ A/ g i are 4U 1626-67 and 
OAO 1657-415. These sources can be monitored to search for 
glitches and anti-glitches with the currently operating X-ray in¬ 
struments and represent good targets for future missions devoted 
to X-ray timing such as the proposed LOFT. 
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